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ABSTRACT. 

Analytic predictions for resonant transport in three generic structures are presented. 
For a structure comprising a normal (N) contact - normal dot (NDOT) - superconducting 
(S) contact, we predict that finite voltage, differential conductance resonances are 
destroyed by the switching on of superconductivity in the S-contact. In the weak 
coupling limit, the surviving resonances have a double-peaked line-shape. Secondly, we 
demonstrate that resonant Andreev interferometers can provide galvonometric magnetic 
fiux detectors, with a sensitivity in excess of the fiux quantum. Finally, for a 
superconducting dot (SDOT) connected to normal contacts (N), we show that the onset 
of superconductivity can increase the sub-gap conductance, in contrast with the usual 
behaviour of a tunnel junction. 
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Recent advances in the fabrication of nanoscale structures have led to increasing 
interest in transport through resonant tunnel junctions and quantum dots [1,2]. In part 
this is due to the new physics associated with Coulomb blockade [3] and in part due to 
growing interest in quantum chaos [4- 7]. The bulk of work in this area has focussed 
on normal structures, but more recently attention has turned to hybrid structures 
involving a superconducting component. It has been demonstrated experimentally that 
the energy gap of a superconducting dot is directly observable through Coulomb blockade 
experiments [8] and theoretical work on incoherent transport through such dots has been 
carried out [9]. In this Letter, we describe the effect of superconductivity on phase- 
coherent transport through resonant structures, in the limit that charging effects can be 
ignored. This limit should be experimentally accessible, because even intimate contact 
with a superconductor will not broaden states below the gap. 

A range of new phenomena and fundamental problems involving coherent transport 
through resonant superconducting hybrids are expected to manifest themselves in a small 
number of generic structures. One such example is a "N-NDOT-S" structure comprising 
one or more normal (N) current carrying leads, in contact with a normal zero dimensional 
"dot" (NDOT), which in turn makes contact with a superconducting (S) lead. In contrast 
with the zero voltage limit, where a general multi-channel description of this structure is 
available [10], there currently exists only a single one-dimensional study of finite voltage, 
resonant transport [11], in which a 5-function potential well, with a single localized state 
is introduced into a one-dimensional insulating barrier. A key prediction of this Letter is 
that finite- voltage conductance resonances are destroyed by the onset of superconductivity 
in such a structure. This is illustrated in figures 1 and 2, which show for a two-dimensional 
system, the conductance of a resonant structure plotted against the Fermi energy of the 
dot. Figures la and 2a show the conductance of a N-NDOT-N structure at zero and finite 
voltages respectively, while figures lb and 2b show the conductance of the corresponding 
N-NDOT-S structure. 
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A second class of structures involves two (or more) separate superconductors S and 
S', with respective order parameter phases (/>'. The resonant energies of a N-SS' or 
N-NDOT-SS' composite will vary periodically with the phase difference (j) — 0' [15-26] and 
in what follows we describe the resonant properties of such structures. A third class is 
formed when a normal lead makes contact with a superconducting dot (SDOT), which in 
turn makes contact with a second normal lead. In such structures, we demonstrate that 
the sub-gap conductance can increase when superconductivity is switched on, in contrast 
with the conventional behaviour of a N-S tunnel junction. 

In what follows, we adopt a multi-channel scattering approach, based on current- 
voltage relations for phase-coherent scatterers written down in reference [27] , which have 
been extended and re-derived in several papers [28-30]. In the absence of inelastic 
scattering, dc transport is determined by the quantum mechanical scattering matrix 
s{E,H), which yields scattering properties of quasi-particles of energy E, incident on 
a phase-coherent structure described by a Hamiltonian H. If the structure is connected 
to external reservoirs by open scattering channels labelled by quantum numbers n, then 
this has matrix elements of the form Sn,n'{E., H). The squared modulus of Sn,n'{E., H) 
is the outgoing flux of quasi-particles along channel n, arising from a unit incident flux 
along channel n' . For channels belonging to current-carrying leads, with quasi-particles 
labelled by a discrete quantum number a (a = -|-1 for particles, —1 for holes), it is 
convenient to write n = (i, a, a), where i labels the leads and a labels all other quantum 
numbers associated with the channels. Then as shown in [27,28], transport properties are 
determined by the quantity 

a,b 

which is referred to as either a reflection probability (z = j) or a transmission probability 
{i ^ j) from quasi-particles of type (3 in lead j to quasi-particles of type a in lead i. For 
a ^ f3, P°^f{E, H) is referred to as an Andreev scattering probability, while for a = 
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it is a normal scattering probability. 

The starting point for our description of resonant Andreev scattering is the following 
general formula for the transition probability between two different scattering channels 
of an open vector space A representing the leads, attached to a closed sub-space B 
representing the scatterer. For n ^ n', the result, which we derive[31] for normal leads 
described by a real Hamiltonian, is 



Tnn' = 4 Trace 

where 



T{n)GBBr{n')GBB^ (1), 



Gbb~^ =gB~^ -(t' -a + iT (2) 

and the trace is over all internal levels of S. In these expressions, T{n) is a Hermitian 
matrix of inverse lifetimes, T = X]„r(n), a and a' are Hermitian self-energy matrices 
and (7b is the retarded Green's function of sub-space B when Hi = 0. The above result 
has been cast in a form which resembles the Breit-Wigner formula [12,13], but is very 
general and makes no assumptions about the presence or otherwise of resonances. 

During the past decade, the Breit-Wigner formula has been applied to a variety of 
problems involving resonant transport in normal-state structures [32-34]. For a normal- 
metallic conductor, under resonant conditions, where the level spacing is much greater 
than the broadening, Biittiker has presented a multi-channel derivation of the Breit- 
Wigner formula through a single resonant level [35] . This limit is recovered from equation 
(1) by restricting the trace to a single level. In what follows, we shall encounter situations 
in which, due to particle-hole symmetry, degenerate states can simultaneously resonate 
and therefore the more general formula (1) is required. 

Consider now a N-NDOT-S structure, where the sub-space B describes the NDOT. 
At zero energy, if the isolated normal dot is on-resonance, then particle-hole symmetry 
ensures that a degeneracy occurs. Hence in this example, the trace in equation (1) can be 
restricted to two terms and equation (2) reduces to an expression involving 2x2 matrices. 
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This yields 



where e± are the particle and hole levels closest to the quasi-particle energy E, d = 

{E -e+ +zr++)(E-e_ -S +tr ) - |(J^_P and we have written S = a-' + u. 

In the presence of a single normal lead, this yields for the electrical conductance in units 

oi2e^/h [14,27,28], 

G = 2P-+(E,iy) = ^£±±£^^^ (4) 

First consider the zero energy limit (E=0), where particle- hole symmetry implies 
that the two levels closest to = satisfy e_ = — e+ and therefore a vanishing particle- 
level e+ is accompanied by a degeneracy. In this limit, S = — S++ and F = F_|__|_. 

Hence 

C - (r. 
((e+ + S++)2+F2.^ + K_P)2 

which was obtained in reference [10] for resonant transport at zero energy. If |cr^_|^ = 
F^_i_, then a resonance will occur when e+ + E+_|_ = 0. Since this involves only a single 
condition on the e+, one expects resonances to occur with approximately the same 
probability when the S-contact is replaced by a N-contact. 

At finite energies, this result is drastically modified, because a resonance can now 

occur only if both {E — e+ — = 0) and — e_ — S = 0). The probability of 

simultaneously satisfying both of these conditions is small and therefore we predict that 
the breaking of the particle-hole symmetry at £^ 7^ destroys finite-voltage conductance 
resonances. From the form of the quartic energy denominator \d\'^, the small number of 
surviving resonances will have a non-Lorentzian line-shape. This is illustrated in figure 
3, which shows plots of the Andreev refiection coefficient Ra = P^^{E,H), for various 
values of coupling to the normal lead; figure 3(a) has the strongest coupling and 3(d) 
the weakest. All quantities are plotted as functions of quasi-particle energy E. Hence 
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as well as a drastic reduction in the probability of finding finite-voltage resonances, we 
predict that when a resonance does occur, the usual Lorentzian line-shape is replaced 
by a double peaked structure, with relative peak heights determined by the difference 
between r+ and r_ at finite energies. The above predictions are confirmed by the results 
shown in figures 1 and 2, which were obtained from an exact numerical solution of the 
Bogoliubov - de Gennes equation, as outlined in reference 28. The inset of figure 2b shows 
numerical results for the fine scale structure of a surviving resonance, in agreement with 
the analytic prediction of figure 3. 

When a composite N-SS" or N-ND0T-S5" structure is formed from two or more 
superconductors, with different order parameter phases, or from a single superconductor 
with an imposed phase gradient [36], transport properties can be significantly modified 
if the phase difference between two points is varied by 27r[15-22]. In experimental 
realizations of such structures[23-26], the phase difference between two superconducting 
contacts is modulated by connecting the superconductors to a macroscopic, external 
superconducting loop, whose phase is controlled by an applied magnetic field. Sub-gap 
quasi-particles can penetrate only a distance of order the superconducting coherence 
length into the superconductor and therefore apart from controlling the phase, the 
macroscopic loop plays no role in determining the s-matrix of the region near the contacts. 
Since the electrical conductance is a periodic function of the phase difference rj, with 
period 27r and since i] changes by 27r when the flux $ through the macroscopic control- 
loop changes by a flux quantum ^q, such structures are galvanometric detectors of flux, 
with a sensitivity comparable with that of a SQUID. We now highlight generic properties 
of resonant interferometers and predict that flux sensitivity can be significantly enhanced. 

The starting point for this analysis is a sub-space B containing one or more 
superconductors, with eigenstates \fi,) and eigenvalues e^, which are periodic functions 
of some dimensionless parameter rj, with period 2tt. In this example, since particle-hole 
degeneracy is lifted by intimate contact with the superconductors, the trace in equation 
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(1) reduces to a single term and the electrical conductance takes the form 

r(n) = 8T,iri+)T,iri-) 

Consider now the situation in which, at = 770, the resonance condition E — Sjy(?7o) = 
is satisfied, where f^viji) = ^v^Vi) + '^vijl)- Then expanding equation (6) about ryo yields 

^/ X ^ 8r^(?7o+)r^(?7o-) 

[aE.(^o)/ar7o]2[r7-r7o]2 + r2(r7o) ^ ^' 

This demonstrates that with varying -q, G exhibits a Lorentzian resonance of width 

r4r/o)/[aS,(r7o)/5r/o]. 

For the case ry = 27r$/$07 noting that Tj^ijf) can vary by at most an amount of order 
Aq as T] varies by 27r, yields an upper bound for [dT,^{r}Q) / dr}Q\ of order Ao/27r. Hence 
in terms of the flux through the external control loop, the resonance width is greater 
than or of order 5$ = 2ti qT ^{rjo) / For simplicity in the above analysis, we have 
considered only a single resonance and a normal lead with no closed channels; the latter 
merely shifts the position of the resonance, while the former may lead to the appearance 
of several resonances per flux quantum. If the temperature T is greater than Ti,{r]Q)/kB, 
then the resonance width will be of order 5$ = 27r$o^B7'/Ao. For a device operating 
at 1 Kelvin, formed from a cuprate superconductor with a transition temperature of 100 
Kelvin, this yields 5$ ~ $o/20. 

Finally we consider a superconducting dot with a uniform order parameter, connected 
to two normal leads. The eigenstates of a superconducting dot satisfy Hslfu) = ^v\fv) 
where Hb is the Bogoliubov-de Gennes operator for the isolated dot. For a dot 
with a uniform real order parameter Aq, if |0) is an eigenstate of the normal dot 
satisfying -ffol^) = e°|</') then the solutions of the Bogoliubov equation are of the form 
6. = ^(i«FTA2, \U) = = (^^^j^l^, where \ut\' = (1 + 6°/6.)/2 and 

= (1 — e^/ejy)/2. At zero energy, this yields for the two-probe conductance derived 
by Lambert [3 1,1 8], 

G = — — (8), 

|-6, + a|^ + z(ri + r2)|2 
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where is the broadening due to contact with lead i. 

If 5e is the level spacing of the normal-state dot, then equation (8) reveals that for 
I Ao| > 5e^ (Fi + F2) the contribution to the conductance from a single level is of order 
G = Hence for a large enough value of Ao/5e, all resonances will be suppressed 

and switching on superconductivity will typically decrease G, a behaviour which is well- 
known for N-S tunnel junctions. However from the form of the denominator in equation 
(8), it is clear that for small values of Aq, the switching-on of superconductivity can 
produce an anomalous increase in the conductance, provided a > e'^. Since will 
be randomly spread between — ^ and ^ this suggests that the probability of a positive 
change is approximately , a result which we have confirmed through numerical solution 
of the Bogoliubov equation for such structures [31]. 

We have presented a theoretical framework and general formulae for resonant 
transport through hybrid normal-superconducting nanostructures. For N-NDOT-S 
structures, we predict that finite- voltage resonances will be almost completely suppressed 
by the switching on of superconductivity and those that survive can have a double-peaked 
line-shape. Such non-Lorentzian resonances have been discussed in other contexts [39] and 
may generate non-exponential delay-time curves. This destruction of resonances implies 
that the ensemble averaged conductance decreases with increasing bias, a behaviour 
reminiscent of zero-bias anomalies in the sub-gap conductance of superconducting- 
semiconducting junctions[40,41]. We have also demonstrated that N-SS' structures can 
possess Lorentzian resonances on a scale much smaller than a flux quantum, which 
suggests that these may provide a new class of magnetometers with a sensitivity at 
least matching that of present-day SQUIDs. Finally we have shown that switching on 
superconductivity in N-SDOT-N structures can produce anomalous positive changes in 
the conductance. 
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Figure Captions. 



Figure 1. For quasiparticles of energy E = Q, the top graph shows the conductance G when 
Aq = 0, as a function of the mean diagonal element eo of a 2-dimension tight binding 
Hamiltonian describing the NDOT. The band width of this Hamiltonian is 8 and the 
Fermi energy is 4 — eo- The lower graph shows corresponding results for Andreev 
reflection coefficient Ra, when the order parameter assumes a non-zero value. 

Figure 2. As for figure 1, except that the energy now takes a non-zero, sub-gap value. The 
inset of figure 2b shows the fine structure of a typical surviving resonance. 

Figure 3. Figures (a) to (d) show plots of equation (3) against energy E, for decreasing values 
of the coupling to the normal lead. 
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